Chaotic mean wind in turbulent thermal convection and long-term correlations in 

solar activity 
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It is shown that correlation function of the mean wind velocity in a turbulent thermal convection 
(Rayleigh number Ra ~ 10^^) exhibits exponential decay with a very long correlation time, while 
corresponding largest Lyapunov exponent is certainly positive. These results together with the 
reconstructed phase portrait indicate presence of a chaotic component in the examined mean wind. 
Telegraph approximation is also used to study relative contribution of the chaotic and stochastic 
components to the mean wind fluctuations and an equilibrium between these components has been 
studied. Since solar activity is based on the thermal convection processes, it is reasoned that the 
observed solar activity long-term correlations can be an imprint of the mean wind chaotic properties. 
In particular, correlation function of the daily sunspots number exhibits exponential decay with a 
very long correlation time and corresponding largest Lyapunov exponent is certainly positive, also 
relative contribution of the chaotic and stochastic components follows the same pattern as for the 
convection mean wind. 

PACS numbers: 96.60.Q-, 47.55.pb, 47.52.+j 
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INTRODUCTION 

Turbulent thermal convection is active background for 
any solar activity. Therefore, the large-scale (coherent) 
circulations, also known as mean winds, spontaneously 
generated by turbulent thermal convection must play a 
significant role in large-scale solar activity. In recent 
decade a vigorous investigation of statistical properties 
of the thermal mean winds has been launched [l| , [3| (see 
for a recent review Q). The winds' dynamics turned 
out to be very complex and many surprising features 
were discovered in laboratory experiments and numer- 
ical simulations. Simple stochastic models were 
replaced by more sophisticated stochastic models 
which also addressed significant three-dimensional nature 
of the phenomenon. In these models, interaction between 
the large-scale wind and the small-scale turbulence pro- 
vides a phenomcnological stochastic driving term. Also 
some deterministic models with chaotic solutions were 
recently suggested QiIIIj and the idea that a large-scale 
instability in the developed turbulent convection (caused 
by a redistribution of the turbulent heat flux) can be an 
origin of the large-scale coherent structures received cer- 
tain experimental support (see Ref. 9] and references 
therein) . 

It is a difficult problem to distinguish between stochas- 
tic and chaotic processes, especially having a developed 
turbulence as a small-scale background. Observation of 
the Lorentzian spectra for the mean winds (reported in 
the Ref. Q) makes this problem even more difficult. It 
is shown recently (see, for instance, Ref. [13] ) that many 
properties of the chaotic processes with exponentially 
decaying correlation function (i.e. with the Lorentzian 
spectra) can be reproduced by some classical models of 
stochastic processes (and vise versa). Therefore, for such 
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FIG. 1: Autocorrelation function of a mean wind velocity for 
a Rayleigh-Benard convection laboratory experiment [^.[TTI]. 
The dashed straight line indicates the exponential decay Eq. 
(1). The insert shows a small-time-scales part of the correla- 
tion function defect in In-ln scales. The straight line indicates 
the Kolmogorov's '2/3' power law for structure function. 



chaotic systems stochastic phenomcnological models can 
be rather useful (and relevant). However, it is necessary 
to reveal an underlying chaotic nature of the phenomenon 
in order to understand real cause of the long-term corre- 
lations in such systems. 



CHAOTIC MEAN WIND 

Figure 1 shows a correlation function of a mean wind 
velocity, v{t), measured in a typical thermal convection 
laboratory experiment in a closed cylindrical container of 



2 



0.007 r 

0.006 
0.005 




0.001 - 

I ■ ' • ' ■ ' • ' • ' 

5000 10000 15000 20000 2500( 

t 

FIG. 2: The pertaining average maximal Lyapunov exponent 
at the pertaining time, calculated for the same data as those 
used for calculation of the correlation function (Fig. 1). The 
dashed straight line indicates convergence to a positive value. 
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FIG. 3: Phase portrait reconstructed from noise-reduced time 
series. 

aspect ratio 1 for Rayleigh number Ra ^ 10^^ (see Ref. 

for description of the experiment and of the other 
properties of the mean wind) . The dashed straight line is 
drawn in the figure in order to indicate (in the semi-log 
scales) the exponential decay 

C(t) ~exp-(r/To) (1) 

The correlation time tq ~ 600s is very large in compar- 
ison with the mean wind circulation period Tc ~ 30s. 
The circulation period (turnover time) provides a time 
scale relevant to the phenomenon physics. The correla- 
tion time normalized by this time scale tq/Tc ~ 20. 
The insert in Fig. 1 shows a small-time-scale part of 



the correlation function defect. The In-ln scales have 
been used in this figure in order to show a power law 
(the dashed straight line) for structure function: {{v{t + 
t) - w(t))2) for T < Tc. 

l-C(r) cx ((i;(i + r) -i;(t))2) cxr2/3 (2) 

This power law: '2/3', for structure function (by virtue 
of the Taylor hypothesis transforming the time scaling 
into the space one 0) is known for fully developed tur- 
bulence as Kolmogorov's power law (cf Refs. P.fl^). 

In order to determine the presence of a deterministic 
chaos in the time series corresponding to the velocity of 
the mean wind for the time scales Tc < t < tq, we calcu- 
lated the largest Lyapunov exponent : Xmax- A strong 
indicator for the presence of chaos in the examined time 
series is condition Xmax > 0. If this is the case, then we 
have so-called exponential instability. Namely, two arbi- 
trary close trajectories of the system will diverge apart 
exponentially, that is the hallmark of chaos. At present 
time there is no theory relating Xmax to Tq^ for chaotic 
systems. There is only a tentative suggestion that their 
values should be of the same order. To calculate Xmax 
we used a direct algorithm developed by Wolf et al [l3|. 
Figure 2 shows the pertaining average maximal Lyapunov 
exponent at the pertaining time, calculated for the same 
data as those used for calculation of the correlation func- 
tion (Fig. 1). The largest Lyapunov exponent converges 
very well to a positive value Xmax — 0.0023s~^. 

The ambivalent nature of the turbulence-induced co- 
herent dynamics one can also see in Figure 3. This fig- 
ure shows a phase portrait reconstructed from the noise- 
reduced time series for the mean wind velocity. 

Thus, we can conclude that the thermal wind under 
study exhibit chaotic features and the long-term expo- 
nential correlation (Fig. 1) is presumably related to this 
chaotic behavior. 

STOCHASTIC-CHAOTIC EQUILIBRIUM 

It is shown in Ref. fl3| that simple telegraph approx- 
imation of stochastic signals can reproduce main statis- 
tical properties of these signals. The telegraph approxi- 
mation of signal v{t) can be constructed as following: 

uit) = vit)/\vit)\ (3) 

From the definition the telegraph approximation can take 
only two values: 1 and -1. Figure 4 shows a compar- 
ison between correlation function of the full signal for 
the mean wind velocity v{t) (cf Fig. 1) and correlation 
function of its telegraph approximation u(t). One can 
see very good correspondence between these correlation 
functions. Therefore, certain main statistical properties 
of the full signal can be studied using the telegraph ap- 
proximation in this case as well. 
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FIG. 4: Autocorrelation functions of a mean wind velocity for 
a Rayleigh-Benard convection laboratory experiment 0| (solid 
curve) and of its telegraph approximation (dashed curve). 
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FIG. 5: Cumulative probability P{t) (Eq. (5)) vs. \n{t/Tc). 
The straight line indicates correspondence to Eq. (6). The 
insert shows a linear dependence of P{t) on t/Tc for t < Tc- 



Very significant characteristic of the telegraph signal 
is duration, r, of the continuous intervals (boxes) where 
the signal takes value 1 or value -1. Probability density 
function of the duration (or life) times r: p(t), for the 
telegraph approximation of the mean wind velocity signal 
was studied in detail in Ref. 0. It has a peak at r ~ Tc- 
In the the range Tc < t < tq the probability density 
exhibits a power law 

p(r) ~ c r" (4) 

with the value of the exponent a close to —1. For r > 
To the probability density p{t) decays exponentially, as 
it should be for the random telegraph signal. Figure 5 



shows cumulative probability 

Pit) = / piT)dT (5) 

Jo 

In the semi-logarithmical scales the straight line indicates 
just the power law for p{t) (Eq. (4) with a ~ —1): 

Pit) » P(r,) + c / r-irfr » P(r,) + c Init/Tc) (6) 

where P{Tc) is a constant part provided by the turbulent 
component to the cumulative probability P{t) for tq > 
t > Tc. On the whole, 

P(oo) ~ AP + c ln(ro/Te) (8) 

where 

AP - P{Tc) + / pir) dT (9) 

J To 

i.e. AP is contribution of the stochastic components: 
turbulence and the large-scale random telegraph signal, 
to the total cumulative probability. In our case there is a 
probabilistic equipartition of the stochastic and chaotic 
components to the total cumulative probability P(oo) = 
1: AP ~ 0.5 ± 0.04. If this probabilistic equipartition 
is universal and the constant c in the power law Eq. (4) 
is also universal (at least asymptotically), then 

ln(To/Pc) = « 0.5/c (10) 

c 

is universal as well. From the present data we can esti- 
mate: c « 0.16±0.01, and ln(To/Tc) « 3 (see also below). 

This kind of probabilistic equipartition is usually re- 
lated to statistical restoration of a symmetry. For in- 
stance, the spontaneous appearing of the mean wind 
results in a spontaneous breaking of mirror (or parity 
- P) symmetry. Statistical restoration of the parity 
(P-symmetry) means probabilistic equipartition of the 
uit) = 1 and uit) = -1 events (Eq. (3)). This prob- 
abilistic equipartition indeed takes place in present case 
and it is clear evidence of the statistical restoration of the 
P-symmetry (the directions of the mean wind rotation: 
clockwise and anticlockwise, are statistically equivalent). 
While the P-symmetry arises from space inversion, the 
T-symmetry arises from time reversal (reversibility). In 
the present case both P- and T-transformations result 
in the same mean wind switching: u — > —u. It can be 
just statistical restoration of the combined PT-symmetry 
(statistical invariance under joint action of parity and 
time reversal) that results in the equal probability of 
the reversible (chaotic) and irreversible (random) switch- 
ings. The statistical restoration of the symmetry in a 
mixed stochastic-chaotic motion indicates an equilibrium 
reached between the stochastic and chaotic components 
(SC-equilibrium) . 
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For the closed orbits in the chaotic systems the Bowen's 
theorem 15| is an analogue of ergodic theorem. This the- 
orem provides a basis for the equality of the average over 
the phase space and the average over large-period or- 
bits. Let us now, in the terms of the Bowen's theorem 
(cf. also [16]), consider a phase volume subregion Afl: 
Oc < Ar2 <C ^ (where ftc is phase volume corresponding 
to the upper turbulent scale Tc) and all the periodic tra- 
jectories passing through it. Different closed orbits can 
be distinguished by their period T. Their distribution 
does not depend significantly on location and boundaries 
of the subregion Afl. Moreover, if we consider the peri- 
odic orbits with periods in the interval (T, T + t) passing 
through the region, then we can use an estimate 
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(11) 



(the coarse graining is defined by the partition of phase 
space by the cells fic)- Let us consider an entropy cor- 
responding to the subregion Ail (cf. Ref. [l3|) and use 
estimate Eq. (11) 



S* = /3 In 



a" 



(3 In 



(12) 



where /3 is a constant (see below, and cf. also Ref. [IT] for 
quantum chaos). The logarithmic growth of entropy can 
imply certain self-similarity of joint probabilistic proper- 
tics of the time series v{t) <-> u{t) (see Eq. (3)). Indeed, 
let us consider probability distribution function, p{v,t), 
for the signal v at the impulses (boxes) of length t <t oi 
its telegraph approximation u. For the stationary signal 
the self-similarity has a standard form [l8| 



(13) 



where / is certain function of the argument v/t^ (the 
statistical restoration of the parity has been taken into 
account here). Let us consider an entropy 



S{t) ~ ~ dv p{v,t) lnp{v,t) 



(14) 



Changing the integration variable: v — > v/t^ and using 
equilibrium condition: S{Tc) — 0, one obtains from Eqs. 
(13) and (14): 



S{t) - /? Ht/Tc) 



(15) 



(cfEq. (12)). 

The insert in Fig. 5 shows a linear dependence of P{t) 
on t/Tc for the turbulent regime that can be related to 
an analytic dependence of P{t) on t/Tc for < t/Tc < 
1 (the linear dependence represents the first two terms 
approximation of the Taylor expansion; the small 'jump' 
at the point t = is presumably related to transition 
from laminar to turbulent motion). In a vicinity of the 
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FIG. 6: Autocorrelation function of the daily sunspots num- 
ber for period 1850-1944y. 



SC-equilibrium P{t) is an analytic function of the entropy 
S{t) = f3\n{t/Tc). Applying again the Taylor expansion 
we obtain 

P{t)=P{Tc) + CoS{t) + ...^P{Tc) + c,(3 Ht/Tc) (16) 

where Co = idP/dS)\s=o (cf. Fig. 5 and Refs. 0,0 
for the transformation t — > In f) . 

It is possible that the above consideration is also ap- 
plicable for the large-scale coherent structures observed 
in other turbulent flows (in turbulent boundary layers, 
for instance [20[). If, for instance, at the SC-equihbrium 
the exponent /3 is a global invariant of motion and its 
value is the same for t < Tc and for t > Tc, then, taking 
into account that the dimensionless entropy S{tq) — 1 
one obtains from Eq. (15): ln(ro/Tc) ~ 1/(3 . For 
the Brownian processes /? — 1/2, whereas for the Kol- 
mogorov's scaling: Eq. (2), one has (3 — 1/3. For 
Kolmogorov's turbulence as a background it results in 
estimate: ln(To/Tc) = 3 (cf. above). If the back- 
ground turbulence produces scaling different from the 
Kolmogorov's one, then value of the /3-exponent can be 
different from 1/3. For the Kraichnan's scaling, for in- 
stance [llj, p — 1/4 and, consequently, ln(ro/Tc) = 4 
(i.e. in the case of the Kraichnan's background turbu- 
lence the correlations related to the large-scale coherent 
structures are even more long-term than those generated 
by the Kolmogorov's background turbulence). 



CHAOTIC SOLAR ACTIVITY 

Let us now look at solar activity, which usually mea- 
sured in the suns pot s number (SSN). In the Babcock- 



Leighton scenario [24]-[25[, for instance, the surface gen- 



erated poloidal magnetic field is carried to the bottom 
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C{x) ~ exp(-T/T„) 
T„=670d 
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FIG. 7: Autocorrelation function of the daily sunspots num- 
ber for period 1964-2008y. 
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FIG. 8; The pertaining average maximal Lyapunov exponent 
at the pertaining time, calculated for the same data as those 
used for calculation of the correlation function shown in Fig. 
6. The dashed straight line indicates convergence to a positive 
value. 



of the convection zone by turbulent diffusion or by the 
meridional circulation. The toroidal magnetic field is 
produced from this poloidal field by differential rotation 
in the bottom shear layer. Destabilization and emer- 
gence of the toroidal fields (in the form of curved tubes) 
due to magnetic buoyancy can be considered as a source 
of pairs of sunspots of opposite polarity. The turbulent 
thermal convection in the convection zone captures the 
magnetic flux tubes and either disperses or pulls them 
trough the surface to become sunspots. This means, that 
the chaotic mean winds of the thermal convection should 
play a significant role in generation and transportation of 
the magnetic field, and consequently in the solar activ- 



ity. Therefore, the chaotic nature of these winds should 
imprint itself in the statistical properties of the solar ac- 
tivity. In particular, the long-term correlation present 
in these winds should result in a long-term correlation 
in the solar activity. Indeed, Figure 6 shows correlation 
function of the daily sunspots number (SSN) for period 
1850-1944y (the SSN data were taken from Next 
to this period two 11-year solar cycles were extraordi- 
nary active (see, for instance, Ref. (i^l) and, therefore, 
we excluded them from the consideration. For the last 
four solar cycles (period 1964-2008y) solar activity was 
again stable (though rather strong) and in figure 7 we 
show the daily SSN correlation function for this period. 
The semi-log scales were used in the Figs. 6 and 7 in 
order to exhibit exponential decay, which corresponds 
to straight line in these scales. It should be noted that 
the maximum entropy method was used in these calcula- 
tions, because it provides an optimal resolution even for 
relatively small data sets. The oscillations in these fig- 
ures correspond to the solar rotation period ~ 27 days. 
This period roughly corresponds to rotation at a latitude 
of 26°, which is consistent with the typical latitude of 
sunspots. The apparent straight lines correspond to the 
exponential decay Eq. (1) (cf Fig. 1). The decay rate 
(or correlation time) takes very close values for these two 
time periods: tq ~ 660d for the period 1850-1944y and 
To ~ 670d for the period 1964-2008y. This correlation 
time is considerably larger than the circulation periods 
inside the convection zone Tc, which are of order of the 
solar rotation period. It should be noted that the corre- 
lation time normalized with the circulation time is close 
to the normalized correlation time of the thermal mean 
wind ln(ro/Tc) ^ 3 (cf previous section). Taking into 
account previous consideration one can assume that this 
long-term correlation in the solar activity could be re- 
lated to the chaotic mean winds in the turbulent solar 
thermal convection. In order to support this assumption 
we show in figure 8 the pertaining average maximal Lya- 
punov exponent calculated for the daily SSN data for the 
period 1850-1944y. The largest Lyapunov exponent con- 
verges very well to a positive value Xmax — 0.0047<i"^, 
indicating chaotic nature of the long-term correlations in 
the solar activity. 



DISCUSSION 

Most of known chaotic systems exhibit exponentially 
decaying correlation functions (see, for instance, Refs. 
[lo|.p8|) or exponentially decaying spectral functions 
(see, for instance, [2^,[33]). It is shown in a recent 
paper Ref. 3l| that solar activity at decadal time scales 
exhibit chaotic behavior of the second kind (see for 
a possible underlying physical mechanism Ref. [s^). 
In present paper we have provided certain evidences 
that at daily (to annual) time scales the solar activity 
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exhibits a chaotic behavior of the first kind. In this 
case the underlying physical mechanism can be related 
to chaotic nature of the mean (coherent) winds sponta- 
neously generated by turbulent thermal convection. The 
observed long-term exponentially decaying correlation 
functions and the positive largest Lyapunov exponents 
for both: mean wind and solar activity, situations can 
be considered as an argument for this suggestion. One 
should also take into account, that the turbulent thermal 
convection plays crucial role in solar activity. Therefore, 
the large-scale (coherent) processes generic to the 
turbulent thermal convection must imprint their statisti- 
cal properties on the solar activity by one way or another. 

The author is grateful to J.J. Niemela, to K.R. Sreeni- 
vasan, and to the SIDC-tcam, the World Data Center for 
the Sunspot Index, the Royal Observatory of Belgium for 
sharing their data and discussions. 
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